Semiclassical approach to quantum Loschmidt echo in deep quantum regions: from 
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In this paper, we study validity of semiclassical predictions for the decay of quantum Loschmidt 
echo in deep quantum regions. This is done numerically in several models (the chaotic regime in 
the sawtooth model and in the kicked rotator model, as well as the critical region of a ID Ising 
chain in transverse field), in which appropriate effective Planck constants can be introduced. In 
particular, we study the change from validity to breakdown of the semiclassical predictions, when 
the effective Planck constants are increased. Our numerical simulations show that the semiclassical 
predictions may work even when the effective Planck constant is not small, in particular, in the 
so-called Fermi-Golden-rule regime. 

PACS numbers: 03.65.Ta; 03.65.-w; 03.65.Yz 



I. INTRODUCTION 

The semiclassical theory is powerful in dealing with 
many problems in various fields of physics [H, ■ It is 
usually expected to work in the semiclassical regime, in 
which the (effective) Planck constant is sufficiently small 
in certain relative sense. An interesting question is how 
deep in the quantum regime could semiclassical predic- 
tions remain valid. Obviously, the answer should depend 
on the physical quantity of interest. 

In this paper, we study a quantity, to which the semi- 
classical approach has recently been found quite success- 
ful in the semiclassical regime. It is the so-called quan- 
tum Loschmidt echo (LE) [H-Q, which is given by the 
overlap of the time evolution of the same initial state 
under two slightly different Hamiltonians, 

M(t) = \m(t)\ 2 with 

m(t) = (* |exp(iff 1 V?i)exp(-iFot/S)l*o), (1) 

where Hi = Hq + eV with e a small quantity and V a 
generic perturbation. The quantity m(t) is usually called 
the amplitude of the LE. The LE gives a measure to the 
stability of quantum motion under small perturbation. 

The LE has quite rich behaviors, depending on the 
nature of the dynamics of the underlying classical sys- 
tem, as well as on the perturbation strength. Usually, 
the LE has a quadratic decay within certain initial time 
interval, as predicted by the first-order perturbation the- 
ory Q. After the initial time interval, in a chaotic sys- 
tem, roughly speaking, the LE has a Gaussian decay 
[El S HI below a perturbative border and has an expo- 
nential decay in the so-called Fcrmi-golden-rulc (FGR) 
regime above the perturbative border with intermediate 
perturbation strength [9j-fl4(. With further increasing 
perturbation strength, in the so-called Lyapunov regime 
under relatively strong perturbation, the LE usually has 



a perturbation- independent decay [TTl fT^ - fl9l ] ; but, in cer- 
tain cases, a perturbation-dependent oscillation in the 
decay rate may also appear [lj, H3 - i20j . 

On the other hand, in regular systems, in the case of 
one degree of freedom, the LE has a Gaussian decay [Ioj |. 
followed by a power-law decay [U, [22[. Meanwhile, in 
the case of many-degrees of freedom, the LE usually has 
an exponential decay for times not very long [23J. 

The above discussed semiclassical predictions for the 
LE have been tested in the deep semiclassical regime, 
with sufficiently small effective Planck constants. Here, 
we are interested in the extent to which they may still be 
valid in the opposite deep-quantum regime, when effec- 
tive Planck constants are not so small. To achieve this 
goal, it is necessary to mainly rely on numerical simu- 
lations in concrete models. Our numerical simulations 
in the sawtooth model and the kicked rotor model show 
that the semiclassical predictions may work well even in 
the deep quantum regime. 

We also study the LE decay in the vicinity of the quan- 
tum phase transition (QPT) of a ID Ising chain in trans- 
verse field. At a QPT, manifested in level crossing of 
the ground level with other level(s), the ground state has 
drastic change(s) in its fundamental properties 24j. The 
properties of QPT have drawn wide attention and have 
been studied extensively in the recent years. In addition 
to traditional quantities employed in the study of phase 
transitions, quantities borrowed from the quantum infor- 
mation field have also been found useful in characteriz- 
ing QPT, such as fidelit y a s the overlap of ground states 
IH-I23] and the LE [H^3li. In the neighborhood of the 
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critical point of the Ising chain, an effective Planck con- 
stant may be introduced and, as shown in Ref . [23| . the 
semiclassical theory may be used to predict the decaying 
behavior of the survival probability, which is a special 
case of the LE. In this paper, we give further investiga- 
tion, in particular, we study validity of the semiclassi- 
cal predictions when the effective Planck constant is not 
small. 

The paper has the following structure. In Sec. II, we 
recall the semiclassical approach to the LE decay. In 
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Sec. Ill, we study the LE decay in the deep-quantum re- 
gion of the sawtooth model and of the kicked rotator 
model. Section IV is denoted to the study of the LE de- 
cay in the vicinity of the QPT of the ID Ising chain in 
transverse filed. Finally, conclusions and discussions are 
given in SecfVl 



II. SEMICLASSICAL APPROACH TO THE LE 

In this section, we recall semiclassical predictions for 
the decay of LE. As well known, the quantum transition 
amplitude from a point ro to a point r in a d-dimensional 
configuration space within a time period t can be ex- 
pressed in terms of Feynman's path integral [32l. f33| . In 
the semiclassical limit, one may use the stationary phase 
approximation to approximately compute the transition 
amplitude. Contributions from paths close to classical 
trajectories give the following well-known semiclassical 
evolution, in terms of Van Vleck-Gutzwiller propagator 
K sc (r;r ;t), 

* sc (r;i)= J dr K sc (r;r ;t)^ (r ), (2) 

where K sc (r; r ; t) = K s (r; r ; t) and 



K s (r; r ;t) 



1/2 



■ exp 



1 



(3) 

Here, the subscript s indicates classical trajectory, 
C\ 12 = \det{d 2 S s /dr l0 d rj )\, (Us is the Maslov index 
counting conjugate points, and S s (r;r ;t) is the action, 
i.e., the time integral of the Lagrangian along the trajec- 
tory s, S s (r; r ;t) = J Q dt'C. 

Let us consider an initial narrow Gaussian wave packet, 



*o(r ) 



1 

— J exp 



:P0 



(ro - ?o) 2 



2e 



(4) 



where (r ,p ) indicates the packet center and £ is the 
dispersion. Semiclassically, the LE is written as M sc (t) = 
\m sc (t)\ 2 , where 



m sc (t) 



dr[*^r;t)]**£°(r;t). 



(5) 



As shown in Refs.fla. l34j|. the amplitude m sc (t) has the 
following explicit expression, 



f \ 



, [' AW - .s (Po-Po) 2 

Po exp [^ A,s ( p °' r ° ; ^ {h/Ci 2 

(6) 

where AS^pg, f ; i) is the action difference along two 
nearby trajectories in two systems H\ and Hq. In the 
first-order classical perturbation theory, with the differ- 
ence between the two trajectories neglected, one has 



A5( Po ,f ;0 -e / dt'V[ Po (t% 
Jo 



(7) 



The LE amplitude in Eq.© can be written as an in- 
tegration over AS, hence, the same for the LE, 



M(t) 



dASe lAs ^ h P(AS) 



(8) 



where P(AS) is the distribution of the action difference 
defined by 



P(AS) ~ ^ J J dp Q S[AS - AS(p Q ,r Q ;t)] 
(Po - Po) 21 



exp 



(9) 



In a chaotic system, averaged over initial states, the 
distribution P(AS) is usually close to a Gaussian dis- 
tribution. When the perturbation is not strong, in the 
so-called FGR regime, AS are relatively small and de- 
viation of P(AS) from the Gaussian distribution is not 
large. In this case, the semiclassical theory predicts the 
following FGR decay for the LE, 



M sc (t) 



-2<r 2 R(E)t 



(10) 



where a = e/h and R(E) is the classical action diffusion 
constant^, 



R{E) 



dt((V[r(t)]-V)(V[r(0)]-V)). (11) 



With increasing perturbation strength, deviation of 
the distribution P(AS) from the Gaussian form can not 
be neglected and one enters into the so-called Lyapunov 
regime. In this regime, when the time is not short, AS 
close to stationary points give the main contribution to 
the LE and the stationary phase approximation predicts 
the following perturbation- independent decay [H[ , 



M sc (t) ocexp[-Ai(t)t], 



where 



Ai(i) 



- — lim In 

t 6x(p)->0 



Sx(t) 



Sx(0) 



(12) 



(13) 



with average taken over initial states. One should note 
that A\(t) is usually not equal to the Lyapunov exponent 
Al, 



Xl = lim — lim In 

t-S-oo t &x(0)^0 



5x(t) 



5x(0) 



(14) 



due to local fluctuations. When the time t is sufficiently 
long, the LE has a decay determined by the long time 
limit of Ai(i) discussed in Ref . (lH| . In a system with a 
homogeneous phase space, i.e., with constant local Lya- 
punov exponent, Ai(t) is given by the Lyapunov expo- 
nent and the LE has the Lyapunov decay [15| . 
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In a ID regular system with periodic classical motion, 
the LE has the following semiclassical expression [22| , 

M sc (t)^c (l + et 2 )- 1/2 e- rt2/{1+et2 \ (15) 
where cq ~ 1 and 



1 / ew p dU \ 

2 V h dpoj 



ew 2 p d 2 U 



2h dp 2 



(16) 



with the derivatives evaluated at the center of the ini- 
tial Gaussian wave packet. Here, w p is the width of the 
initial Gaussian wave packet in the momentum space, 
U = y- J Q " V dt and T p is the period of the classical 
motion. 

On the other hand, in a regular system with many 
degrees of freedom, when the time is not long, due to 
the large number of frequencies, the system has a motion 
similar to a chaotic one (though the motion shows quasi- 
periodicity for sufficiently long times). As a result, the 
LE has a FGR-type decay in Eq. (fT0| 
not very long j23| . with 



when the time is 



R(E) 




V{t)dt 




(17) 



Finally, we note that Eq. © is a general expression, not 
restricted to the semiclassical limit. To show this point, 
one may use Feynman's path-integral formulation. For 
brevity, lets us write Feynman's propagator as 



K F (r , r ; t ) = N ^ exp {iS a (r , r ; t)/H} , 



(18) 



where a indicates possible paths going from ro to r within 
a time interval t and J\f is the normalization coefficient. 
Using this propagator, the exact time evolution of the 
wave function ^(r, t) can be written in a form similar to 
that in Eq.Q, with K sc replaced by Kp. Then, substi- 
tuting the expression obtained into the definition of m(t) 
in Eq.([T]), one obtains 

m(t) = NN' J dror'o ^ exp(iAS e /h)M^W (r' ), 

(19) 

where AS e = S^-S^, 1 . It is seen that the LE amplitude 
m(i) can always be written as an integration over the 
exact action difference AS 1 ,,, with the distribution P(AS) 
defined accordingly, as a result, the LE can always be 
written in the form of Eq.([S]). 
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FIG. 1: LE decay in the FGR regime in the sawtooth model, 
for different values iV of the dimension of the Hilbert space. 
Parameters: K = 2.0 and a = 0.5. The LE (averaged over 
initial states) has two stages of decay. In the first stage, it fol- 
lows the semiclassically-predicted FGR decay (solid straight 



lines), namely, e 



decays with slower rates. 



120 



Eq. (|27|) ; in the second stage, it 
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FIG. 2: Variation of td (empty squares) with N in the saw- 
tooth model, where td is the time at which the second-stage 
decay of the LE in the FGR regime appears (see Fig[T|). For 
comparison, we also plot t n (solid circles), the time at which 
the FGR decay is expected to reach the saturation value 1 /N , 
i.e., t n = (451n(A0)/(o- 2 7r 4 ). Parameters: K = 2.0, a = 0.2. 



Two kicked models 



III. LE DECAY IN DEEP QUANTUM REGION 
OF TWO KICKED SYSTEMS 

For the purpose of studying validity of semiclassical 
predictions in the deep quantum region, it would be con- 
venient to employ models in which effective Planck con- 
stants can be suitably introduced. In such a model, the 
value of effective Planck constant gives a natural measure 
for quantum "deepness" . 



We employ the sawtooth model and the kicked rotator 
model, whose Hamiltonians have the following form, 

1 °° 
H=-p t + V(r)^2S(t-nT), (20) 

where V(r) = —k(r — n) 2 /2 for the sawtooth model and 
V(r) = fccosr for the kicked rotator model. Here, for 
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FIG. 3: Similar to Fig[T] but for a = 3.0 in the Lyapunov 
regime of the sawtooth model. The LE has a decay close to 
the semiclassically-predicted Lyapunov decay, namely, e _Al/t 
(solid lines), for N > 2 6 . 



FIG. 4: Variation of the deviation D in Eq.j29j) with N in 
the Lyapunov regime of the sawtooth model, with parameters 
K = 2.0 and a — 3.0. The value of D keeps small for N > 
N c = 64 and becomes large when N is smaller than N c . 



simplicity in discussion, wc consider their dimensionlcss 
form. Hereafter, we take unit Planck constant, H = 1. 

The classical dynamics in the kicked rotator model gen- 
erates the standard map, 

Pn+i = Pn + Ksin(r n ) (mod 2ir), 

r n +i = r n +p n+ i (mod 2vr), (21) 

where p n = Tp n , K — kT. The classical motion is reg- 
ular for sufficiently small K and is almost chaotic for K 
larger than 6 or so. In the sawtooth model, one has the 
following classical mapping, 

Pn+i = Pn + K(r n - 7r) (mod 2w), 

r n +i = r n +p n+1 (mod 2ir). (22) 

Equation. (|22|l can be written in the matrix form 



Pn+l 
r„+i - 7T 



1 K 
1 K + l 



Pn 



(23) 



The constant matrix in the above equation possesses two 
eigenvalues 1 + (K ± \/(K) 2 + AK) /2. The motion of the 
classical sawtooth model is completely chaotic for K > 0, 
with the Lyapunov exponent 



X L = ln({2 + K + [(2 + Kf - 4] 1 / 2 }/2), 



(24) 



given by the larger eigenvalue of the constant matrix. 

We utilize the method of quantization on torus to get 
the quantum versions of the above two classical systems, 
with periodic boundary condition for the coordinate and 
momentum variables, < r < r m , < p < p m [l|, l35h - 
|37| . For a Hilbcrt space with dimension N, an effective 
Planck constant can be introduced, denoted by h e s, with 
^cff = TmPmjN . In the specific choice of r rn = p rn = 2ir, 
which will be taken in what follows, we have S e ff = 2n/N. 



The value of ft e ff gives a measure for the "deepness" in 
the quantum region, with relatively-large value of h c g 
implying being relatively deep in the quantum region. 
The evolution operator for one period of time T, with 
T = 2tt/N = h c g, is written as 



U = cxp 



~2 



2h e g 



cxp 



— 7(r) 

fteff 



(25) 



where V(r) = TV(r). 

In the two kicked models discussed above, the quantity 
R(E) appearing in the FGR decay in Eq. ([T0f has the 
following expression [H, [39| , 



i 00 



(26) 



where C(l) = ({V[r(l)] - (V)}{V[r(0)] - (V)}). In the 
sawtooth model with an integer K, C(0) = 7r 4 /45 and 
C(l) = for I ± 0, hence, 



M sc (t) 



(27) 



Meanwhile, in the kicked rotator model, R(E) is a func- 
tion of the parameter K and does not have an explicit 
analytical expression. 



B. Numerical results in sawtooth model 

In this subsection, we discuss our numerical simula- 
tions obtained in the sawtooth model. In the FGR 
regime, when N is not large, it was found that, be- 
yond some initial times, the LE has two stages of de- 
cay (see FigJTJ) : In the first stage, the LE follows the 
semiclassically-predicted FGR decay, while, in the sec- 
ond stage, it is somewhat slower than the FGR decay. 
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After these two stages of decay, the LE oscillates around 
its saturation value, which is on average approximately 
equal to 1/N [13] • Interestingly, the first-stage decay of 
the LE exists even for small dimension N of the Hilbcrt 
space, hence, for large value of the effective Planck con- 
stant. 

We use td to indicate the transition time of the above- 
discussed two stages of decay of the LE, i.e., the time at 
which obvious deviation from the FGR decay appears. 
As seen in FigJ5J td increases with increasing N. In ad- 
dition, we observe that the second-stage decay of the LE 
approaches the FGR decay with increasing N, that is, 
the difference between the decay rates in the two stages 
decreases. 

To get further understanding in the above-discussed 
first and second-stage decay of the LE, let us reconsider 
the expression of the LE in Eq.©. The distribution 
P(AS) always has some deviation from its Gaussian ap- 
proximation, which we denote by Pq with G standing for 
Gaussian, that is, 



P(AS) = P G + AP. 



(28) 



The above discussed numerical results imply that the de- 
viation AP is not sufficiently large for times shorter than 
td (beyond some initial times), as a result, the LE still 
follows the FGR decay. However, for times beyond td, 
the deviation can not be neglected. 

As discussed previously, Eq.© is not just a semiclas- 
sical expression, but, is an exact expression, if the distri- 
bution P(AS) is appropriately defined in terms of con- 
tributions from Feynman paths. The deviation AP is 
mainly induced by two factors: One comes from contribu- 
tions not included in the stationary phase approximation, 
which is used when deriving the scmiclassical propagator 
from Fcynman's path integral formulation. Meanwhile, 
the other is related to the fact that the right hand side 
of Eq.© for classical trajectories does not have an exact 
Gaussian form. As discussed above, an obvious second- 
stage non-FGR decay of the LE appears for relatively 
small values of the dimension N, which correspond to 
relatively large values of the effective Planck constant. 
This implies that the above-mentioned first factor plays 
the major role here, i.e., the deviation is mainly induced 
by non-semiclassical contributions. 

Next, we goes to the Lyapunov regime in the saw- 
tooth model. In this regime, we did not observe a two- 
stage decay similar to that discussed above in the FGR 
regime. For large N, as shown in previous work [Tij ]. 
the LE has approximately the semiclassically-predictcd 
Lyapunov decay, with decaying rate given by the Lya- 
punov exponent. (The sawtooth model has a homoge- 
neous phase space.) With decreasing N, as shown in 
FigEl the decay of the LE gradually deviates from the 
Lyapunov decay. 

In order to quantitatively characterize the above dis- 
cussed deviation of the exact LE decay from the semi- 
classical predictions in the Lyapunov regime, we have 
studied the standard deviation of x n = | In M e (t = 




2.0 2.4 2.8 3.2 3.6 4.0 4.4 4.8 5.2 



FIG. 5: Variation of N c with a in the Lyapunov regime of the 
sawtooth model with parameter K = 1.0. 



n) — \nM sc (t = n)\, where M e (t) denote the exact nu- 
merical results. Specifically, we have studied the quantity 
D, 



D = 



1 



M 
n=l 



(29) 



where 



M S 



is the average value of x n . Our 



numerical simulations show that the value of D remains 
small for large N and is not small when N is below some 
value, which we denote by N c (see Fig|4j) . That is, the 
scmiclassical prediction for the LE decay works well for 
N above N c , but, not well for N below N c . 

The value of N c was found dependent on the param- 
eter (j, as seen in FigJS] On average, N c increases with 
increasing a. This dependence may be related to a re- 
quirement used in the derivation of the above-mentioned 
scmiclassical predictions for the LE decay, namely, e be- 
ing small. Indeed, due to the relation e = erft c ff = 2tt<j/N, 
to keep e small, larger value of a corresponds to larger 
value of N. 



C. Numerical results in the kicked rotator model 

In the kicked rotator model, numerically we found be- 
haviors of the LE more or less similar to those in the saw- 
tooth model discussed in the previous subsection. Specif- 
ically, in the FGR regime, when N is not large, we also 
observed a two-stage decay of the LE; but, in this model, 
the second-stage decay is faster than the first-stage FGR 
decay. See Fig|6] for some examples, where the value of 
R(E) in the FGR decay was computed numerically, mak- 
ing use of Ea.pp]). 

The kicked rotator does not have a homogeneous phase 
space, hence, in the Lyapunov regime, the semiclassical 
prediction for the LE decay is not given by the Lya- 
punov exponent of the underlying classical dynamics, 
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t t 



FIG. 6: Same as in Fig[T] for the kicked rotator model with 
parameters K = 11.0, a — 0.3 in the FGR regime. The 
semiclassical prediction of the FGR decay, e~ 2<T w jth 
R(E) = 0.375, is indicated by solid line. 



FIG. 7: Similar to Fig(6] but for the Lyapunov regime with 
parameters K = 15.0 and a = 20.5 in the kicked rotator 
model. The semiclassical prediction, given by Eq. (|12[l . is in- 
dicated by solid lines. 



but, is given by Eq. (fT2|) . As expected, only for large N, 
did we numerically find agreement between predictions 
of Eq.(fT2"j) and the exact LE decay beyond some initial 
times. For example, as shown in FigEJ the agreement 
is good for N = 2 13 and t e [4,8]. On the other hand, 
the agreement is not good for N < 2 9 ; for this value 
of N, the LE approaches its saturation value before the 
semiclassically-prcdictcd decay in Eq. (IT2l can be seen. 



IV. LE DECAY IN THE VICINITY OF A 
QUANTUM PHASE TRANSITION 

As discussed in Ref.[23j], the semiclassical theory can 
be used in predicting the LE decay in the close neighbor- 
hood of those QPTs whose ground levels are infinitely 
degenerate at the critical points; the closer the control- 
ling parameter A is to the critical point A c , the better the 
semiclassical theory might work. This is related to the 
largeness of the density of states near the critical point. 

In this section, we study the LE decay close to a crit- 
ical point of a ID Ising chain in a transverse field. As 
discussed below, an effective Planck constant may be in- 
troduced in the low energy region in this model, which 
is inversely proportional to the number N p of the spins. 
The system undergoes a QPT at the critical point in the 
thermodynamic limit N p — > oo. Here, we are particularly 
interested in validity/breakdown of the semiclassical pre- 
diction of the LE decay in the case of finite N. In this 
model, Hq in the definition of the LE in Eq.([T]) is taken 
as H(\ ) and H\ as H(X). The initial state |^o) is cn0 ~ 
sen as the ground state of H(Xo), thus, the LE here is in 
fact a survival probability. 

In this model of ID Ising chain, the dimensionlcss 



Hamiltonian is written as 

H(X) = -J^(afa? +1 + Xaf). (30) 

The spin-spin interaction intends to force the spins to po- 
larize along the z direction, while the transverse field in- 
tends to polarize them along the x direction. The compe- 
tition of the two interactions results in two critical points, 
A c = ±1, with ferromagnetic phase for — 1 < A < 1 and 
paramagnetic phase for |A| > 1. Without the loss of gen- 
erality, we consider the critical point A c = 1. 

The above Ising Hamiltonian can be diagonalized by 
utilizing the Jordan- Wigner and Bogoliubov transforma- 
tions, giving [H, 0, |H(, 

ff(A) -1/2), (31) 

k 

where b k and b k are fermionic creation and annihilation 
operators, e k is the corresponding single quasi-particle 
energy, 

e k = 2v/l + A 2 -2Acos(fca), (32) 

and k = 2Trm/aN p with m = —M, —M+l, • • • , M. Here, 
a is the lattice spacing and N p = 2 M + 1. 

As discussed in Ref.[23j], in the very neighborhood of 
the critical point with A sufficiently close to A c and for 
sufficiently large N pi the low- lying states can be mapped 
to bosonic modes (labeled by a), with single-particle en- 
ergies e b a ~ n a SE, where n a = 1, 2, . . . and SE = 4tt/N p , 
and, furthermore, the semiclassical theory predicts the 
FGR decay in Eq.([IU]) with R{E) given by Eq.(|rT|). The 
form of the single-particle energy e h a suggests that one 
may introduce an effective Planck constant h e s, 

h cS = 5E = Air/Np. (33) 
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0.0 0.4 0.8 1.2 1.6 2.0 



t 

FIG. 8: LE decay of the Ising chain for different values of 
N p , with parameters Ao = A c — 4 x 10~ 2 , A = A c — 10~ 2 and 
A c = 1.0. For large N p and relatively long times, the LE 
has an exponential decay as predicted by the semiclassical 
theory. (A solid straight line is drawn for guiding eyes.) Note 
also that the LE has a good scaling behavior of InM oc —N p t 
for large N p . For relatively small N p (N p = 25), the LE has 
neither the exponential decay nor the scaling behavior. 
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FIG. 9: Variation of the deviation D with the spin number 
N p in the Ising model, with parameters Ao = A c — 4 x 1CP 2 
and A = A c — 1(T 2 . 

Then, the mode energies are written as e b a = h e su> a , with 
io a ~ n a . 

Numerically, we observed the semiclassically-predicted 
exponential behavior of the LE when N p is relatively 
large (see Fig(5|). But, for sufficiently small N p , obvi- 
ous deviation from the exponential decay appears. We 



have further calculated the deviation D in Eq. ([29|) for 
x t = | In M e (f) — lnM sc (i) (sec Fig|9|), where the slope 
of lnM sc was computed in the large N p limit. It is seen 
in Fig(5]that, with decreasing N p , the value of D first 
remains small, then, begins to increase rapidly around 
N p = 100, indicating breakdown of the semiclassical pre- 
diction. 

Furthermore, as seen in FigJSl when N p is large, the 
LE has a good scaling behavior of InM(i) ~ N p . To 
understand this phenomenon, we note that here the per- 
turbation eV in the definition of the LE takes the form 
of e = (A — A) and 

N p 

V = J^<T?. (34) 

i=l 

Hence, according to Eq. (fTT|) , the quantity R(E) is given 
by the square of the summation of N p terms with mean 
zero, each of which is a time integration of e(af — erf)- As 
discussed in Ref.dH, in the low energy region and near 
the critical point of this model, the system has a clas- 
sical counterpart with many degrees of freedom, hence, 
possesses a classical motion with many frequencies. This 
implies that the time integrations mentioned above can 
usually be regarded as being uncorrclatcd. Then, for 
large N p , the quantity R(E) is approximately propor- 
tional to N p , hence, lnM(i) ~ N p . 



V. CONCLUSIONS 

In this paper, we have studied the change from va- 
lidity to breakdown of semiclassical predictions for the 
LE decay in several models, when the effective Planck 
constants are increased and the systems move from the 
semiclassical region to the deep quantum region. Our nu- 
merical results show that some semiclassical predictions 
for the LE decay work well even in the deep quantum 
region. 

In particular, in the FGR regime of quantum chaotic 
systems, there is always some time interval within which 
the LE follows the FGR decay; but, the length of this 
time interval decreases with increasing effective Planck 
constant. Beyond this time interval, deviation from 
the FGR decay appears, which is induced by non- 
scmiclassical contributions. In the Lyapunov regime, the 
semiclassical prediction remains valid when the effective 
Planck constant is not quite small, but, it becomes in- 
valid for sufficiently small effective Planck constant. 
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